A surprising relation between double exchange and Heisenberg model spectra: 

Application to half-doped manganites 
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The Zener polarons recently found in half-doped manganites are usually seen as mixed valence 
entities ruled by a double exchange Hamiltonian involving only correlated electrons of the metals. 
They can however be considered as ferrimagnetic local units if the holes are localized on the bridging 
oxygen atoms as implicitely suggested by recent mean-field ab initio calculations. In the latter 
case, the physics is ruled by a Heisenberg Hamiltonian involving magnetic oxygen bridges. This 
paper shows that the spectra resulting from the resolution of both models are analytically identical. 
This single resulting model spectrum accurately reproduces the spectrum of Zener polarons in 
Pro.6Cao.4MnC>3 manganite studied by means of explicitely correlated ab initio calculations. Since 
the physics supported by each model are different, the analysis of the exact Hamiltonian ground 
state wave function should a priori enables one to determine the most appropriate model. It will be 
shown that neither the spectrum nor the wavefunction analysis bring any decisive arguments to settle 
the question. Such undecidability would probably be encountered in experimental information. 

PACS numbers: 71.27. +a,75.10.Dg,75.30.Et,75.47.Lx 



Strongly correlated materials have attracted the inter- 
est of chemists and physicists for their fascinating macro- 
scopic properties 0,0, 13 • Adequate treatments of collec- 
tive effects as well as the resolution of the microscopic 
origin of these properties are crucial to the understanding 
of these materials. Both alike require the determination 
of accurate effective Hamiltonians. In doped mangan- 
ites as well as in many oxydes, the various phase tran- 
sitions occuring under doping are not fully understood. 
As well, relevant and quantitative effective Hamiltonians 
able to reproduce the microscopic order in the most com- 
plex phases are rarely available. 

A recent crystal determination of Pro.6Cao.4MnC>3 
manganite pj has revealed the presence of Zener polarons 
i.e. trappings of holes (or electrons) within pairs of Man- 
ganese sites. The Zener polaron was initially assumed to 
be ruled by a double exchange model0, Q giving rise to 
an octet ground state of the dimer. In this description 
each hole is delocalized over two Mn sites resulting in a 
resonance between Mn 3+ 2 ~Mn 4+ and Mn 4+ 2 ~Mn 3+ . 
However, recent ab initio calculations on the closely re- 
lated Lao.5Cao.5MnO30,0 show an important O to Mn 
charge transfer, resulting in a localization of the holes on 
the bridging oxygens. This charge distribution suggests 
a dominant purely magnetic local order (Mn 3+ 0~Mn 3+ ) 
in which the Zener polarons would be ferrimagnetic enti- 
ties involving an antiferromagnetic coupling between the 
magnetic oxygen and the high spin manganese centers. 
In the latter case the model Hamiltonian which provides 
a relevant description of the local electronic order is a 
Heisenberg Hamiltonian. The holes localization has been 
a matter of debate in the literature and concerns a large 
number of materials like doped cuprates 8] (where there 
is at most one unpaired electron per metal resulting in 
a t-J model), nickelates or manganites (in which there 



are several open-shells per center leading to a double 
exchange interaction). Besides, as far as the correlated 
electrons are only localized on the metallic centers, the 
comparison of the models has been the subject of in- 
tense discussion and controversy which ended up with the 
demonstration that the two models were strictly different 
and lead to different spectra^. The aim of this paper 
is to show that when considering a magnetic oxygen in 
the Heisenberg model the analytical resolution of the two 
Hamiltonians leads to identical spectra. For this purpose 
we have derived an analytical solution of the Heisenberg 
Hamiltonian energies for three magnetic centers. 

Since the discrimination between the two models is im- 
possible from the spectrum determination, we will then 
adress this question from wavefunction analysis. Indeed, 
the projection of the ab initio ground state wavefunc- 
tion onto both model spaces should a priori determine 
the choice. However, it will be shown that when the 
charge of the oxygen is around -1.5, the wavefunction 
analysis does not infer any decidable argument. In or- 
der to illustrate this assertion, we have extracted the two 
corresponding models from the correlated ab initio calcu- 
lated spectrum of the Pro.6Cao.4Mn03 manganite Zener 
polaron and analyzed both its spectrum and its ground 
state wave function. 

The here-developped solutions are analytical and could 
be applied to any material having a bridging ligand 
and an odd number of correlated electrons per metallic 
dimer. Thus, the here-developed equations have been 
generalized to any number n of open-shells per cen- 
ter. The first considered model, is the Zener double ex- 
change model^(J including the Heisenberg-type antifer- 
romagnetic contribution proposed by Papaefthymiou and 
Girerd|9llTl|. It is likely to describe configurations com- 
patible with a M(d( n - 1 / 2 ))0 2 -M(d n - 1 / 2 ) mixed valence 
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Heisenberg model space is also built from products of 
atomic ground states. The S 2 operator can be written 
as 
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S 2 = (S, + S,. + S„V (2) 

2(S;.So + S r .So + S r .S;) 

Replacing S r .Si by its expression as function of S 2 in 
the Hamiltonian, one gets 
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where S' = S r + S;. As S'.S = |(S 2 - S' 2 - S 2 ) 
Hamiltonian can be rewritten as: 
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FIG. 1: Sample cases of electronic distributions (left: double 
exchange model; right: Heisenberg model) illustrating the two 
kinds of systems. In a) — nf (Nickelates) while in b) 
n„ — n u + 2 (Manganites) . 



electronic structure resulting from a resonance between 
M(d("))0 2 -M(d™- 1 ) and M(d("- 1 ))0 2 -M(d n ) where n 
and n-1 are the numbers of open shells on the two metal- 
lic centers. The model is based on the idea that the 
spectrum can be reproduced by considering the metal- 
lic ions in their atomic ground states and two effective 
interactions, namely the hopping integral / of the extra 
electron between the two metal-centered orbitals and an 
overall Heisenberg exchange integral J between the elec- 
trons of the other open shells. The derealization of the 
extra electron favors the highest spin state (of total spin 
^max) while the exchange integral stabilizes the low and 
intermediate spin states for an overall antiferromagnetic 
contribution (J < 0). In the particular case of a symmet- 
ric homonuclear bimetallic complex, the eigenenergies of 
the double exchange model are analytically known. This 
Hamiltonian generates two eigenstates for each value of 
the total spin S, the energies of which E Z (S, ±) are given 
by the expression: 



where the zero of energy is taken as the mean value of 
the highest spin states energies. 

If the electronic structure is dominated by 
M(d™)0~M(d™) charge tranfer configurations, 
a Heisenberg Hamiltonian is expected to rule 
the physics of the system. Its specific form is 
H H = — JiS;.So — JiS r .So — ^2S/.S r where Ji is 
the overall spin exchange interaction between the left 
(S;) and right (S r ) metal spin, and J\ the magnetic 
interaction between the metal and the oxygen ions. The 
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Since the oxygen has only one unpaired electron, i.e. 
So = 1/2, the values of S' can only be S' = S± i and an 
analytical expression of the eigenenergies can be found: 
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where S^ ax = Si + S r + 5. For the highest spin multiplic- 
ity state, only the E (5, — ) root has a physical meaning 
since 5' has the single value S' = S — | . 

In order to compare the two models, the same zero of 
energy must be used in equations (IJ and (JSJ . Let us use 
the double exchange model zero of energy. Its expression 
as a function of the Heisenberg model is: 
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The relation between S^ ax and S^ ax is fixed by the 
number n u of unpaired electrons in the two models (n^f 
n z for the Heisenberg and double exchange models, 

= Sg ax . On the 
SLx + ±- These 



and 

respectively). If n% = nf , then S% a 
contrary, if <nfi = nf +2 then S% ax = 
two situations are illustrated in figure ^ The first one 
(schemes a) corresponds, for instance, to dopped nicke- 
lates. The delocalized electron of the double exchange 
model (left scheme) jumps from a doubly occupied or- 
bital towards a singly occupied one. The second situa- 
tion (schemes b) is representative of half-doped mangan- 
ites. The hopping electron of the double exchange model 
jumps from a singly occupied orbital towards the empty 
orbital of the other manganese. 

The Heisenberg energies shifted by the mean value Eq 
of the S z nax states become: 
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Thus, from comparisons between equations and J7J, 
an important conclusion is that the energies of the two 
models are exactly the same if one replaces t and J by 



t = \{Ji-J2)[si + \ 



and J = J-2 



(8) 



Let us notice that when S^ lax 



= Sfnax + 1) the en- 
ergy of the S^ lax state cannot be obtained in the double 
exchange model, since this state does not belong to the 
model space. Its energy is only given by the Heisen- 
berg expression. That means, in dopped manganites 
case, that the double exchange model contains two dou- 
blet, two quartet, two sextet and two octet (S = 7/2) 
states whereas the Heisenberg model also describes a de- 
cuplet state (S = 9/2). On the contrary, when S^ ax = 
Smaw the E H (S^ naxl +) solution is meaningless for the 
Heisenberg model (see eq.|5]and its explanation) whereas 
E (Sm ax , +) is a solution of the double exchange model. 
Nevertheless, we still have E H (S% ax , +) = E z (S^ ax , +). 
So, the complete description of the spectrum is possi- 
ble from the Heisenberg expression while the Heisenberg 
model space does not contain the E H (S^ ax ,+) state. 
One therefore has the unexpected situation in which two 
model Hamiltonians working on different model spaces 
generate spectra which coincide except for one state. If 
Smax = Smax + 1 the double exchange spectrum is in- 
cluded in the Heisenberg one while it is the Heisenberg 
spectrum which is fully described by the double exchange 
model when S% ax = S z ax . 

In order to check the relevance of the two considered 
models on a real material, let us turn now to the ex- 
plicitely correlated ab initio study of the Pro.6Cao.4Mn03 
manganite. Since the choice between the two models 
turns out to be impossible from the only determina- 
tion of the spectrum, the analysis of the ground state 
wavefunction becomes compulsory in order to discrimi- 
nate between them. The M^O^^ cluster (constituted 
of the dimer of Mn sites surrounded by the oxygen 
atoms of their coordination sphere) is analyzed using 
the embedded cluster method. In order to reproduce 
the crystal environment, the cluster is embedded in a 
bath that reproduces its main effects (Madelung and 
Pauli exclusion) 0, 0] . Then the cluster is studied us- 
ing correlated ab initio calculations and large basis sets 
(l3| . MS-CASPT2 calculations Q (the many-body sec- 
ond order perturbation scheme of the MOLCAS codefl5| 



which has been successfully applied to many spectro- 
scopic studies (ljjl are performed (for a more detailed 
description, see references and Gj). The so-obtained 
spectrum constituted of height states compares well with 
the model spectrum 12] . Table U reports the optimized 
parameters and the mean error per states of the two 
models. This error is calculated as the energy differ- 
ence between the computed (ab initio) states and those 
reproduced by the model divided by the spectrum width 
and the total number of states. These results show that 
both models appear to be appropriate to the description 
of the polaron physics (it is interesting to note that a 
dramatic improvment of the model resulting in a 0.7 % 
of error is obtained when the Non-Hund states are ex- 
plicitely considered, as shown in ref. 0])- The large 
value of J\ reveals a partially covalent interaction be- 
tween the Mn atoms and the bridging oxygen. However, 
the energetic position of the decuplct ( 25+1=10) state 
in the narrow vicinity of the highest octet state is in 
perfec t ag reement with the prediction of the Heisenberg 
model|l2j. This result attests of the partially magnetic 
character of the oxygen. This intermediate character of 
the metal-oxygen interaction is confirmed by the shape of 
the symmetry-adapted natural orbitals which are repre- 
sented in figure Et- The derealization between the met- 
als and the oxygen ligand is large in both the essentially 
metallic orbital (1 and 2) and the essentially 2p CT oxygen 
centered orbital (3). The occupation number (1.6) of the 
essentially 2p CT oxygen centered orbital (3) is in between 
the two limit values of 2.0 in the DE model and 1.0 in 
the Heisenberg one. 





12) 




(3) 




(2) 



Ji J 2 t J error (%) 
Heisenberg 0.59 0.07 - - 3.3 
Double exchange - - 1.05 0.07 3.3 

TABLE I: Effective interactions (in eV) of the two models 
extracted from the ab initio spectrum and errors of the corre- 
sponding predicted spectrum compared to the ab initio one. 



FIG. 2: Canonical delocalized orbitals (a); strongly localized 
orbitals (b) and left or right semi- localized orbitals (c). 



In order to discriminate between the two hamiltoni- 
ans we have compared the weight of the ground state 
model wavefunctions in the exact Hamiltonian ground 
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state wavefunction. Model Hamiltonians assume local- 
ized orbitals sets which must be properly defined. The 
symmetry-adapted orbitals set is not appropriate to the 
analysis of the Heisenberg wave functions due to the large 
derealization of orbitals (1) and (3) (see fig. CJO- In order 
to calculate the weight of the Mn 3+ 0~Mn 3+ electronic 
configuration, orbitals (1), (2) and (3) have been relo- 
calized using Boys localizations |l7j . The three strongly 
rclocalized orbitals are represented in Figure 2b. 

The Heisenberg Octet ground state wave function is : 

W^-) = % (Qi) i D + i/ 2 (Q 2 ) I ,-4(Q 2 ) ; D_ 1/2 (Q 2 ) r 

+ (Q 2 ),D +1/2 (Q 1 ) r ] 

where (Q 2 )/D_!/ 2 (Q 2 ) is the product of the local left 
and right Sz = 2 quintet components (Q 2 ) on the two Mn 
sites and a local Sz = —1/2 doublet component D_i/2 on 
the oxygen. (Qi). r D +1/2 (Q 2 ) ; (resp. (Q 2 ) r D +1/2 (Qi) ; ) is 
the product of the right (Qi) r (resp. left (Qi) ; ) quintet 
Sz = 1 component on Mn r (resp. Mn;) and D +1/ / 2 is the 
Sz = +1/2 doublet component on O. 

The projection of the calculated ground state wave 
function onto the Heisenberg ground state gives a weight 
of 0.74 showing the adequation of this model Hamiltonian 
to reproduce the polaron physics. 

The same wave function analysis has been performed 
on the double exchange model space. The double ex- 
change octet ground state wave function is : 

* z (l~) = ^[(Q 2 ) i (Q3/ 2 ) r + (Q3/ 2 ) / (Qi) r ] 

where Q3/2 is the S z = 3/2 component of a local quar- 
tet. This model assumes a two-electron occupied orbital 
centered on the oxygen which is actually well described 
within the symmetry-adapted delocalized picture. The 
double exchange eigenvectors are linear combinations of 
two resonant forms where the extra electron is either on a 
left or on a right localized orbital. The adequate metallic 
orbitals are semi localized on the left and on the right 
metals with correct derealization tails of the bridging 
oxygen. The so-obtained three orbitals are represented 
in Figure 2c. The weight of the calculated wave function 
on the double exchange model space in this orbitals set is 
0.72 (to be compared to 0.74 for the Heisenberg model), 
revealing that both models are equivalently suitable for 
describing the physics of the polaron. 

As a conclusion, since the spectra of the double ex- 
change and the Heisenberg models are identical, the 
only criterium to choose the most appropriate model 
for a specific system should be the wavefunction pro- 
jection onto both model spaces. In limit cases where 



the oxygen charge is close to -1 (or -2), i.e. for domi- 
nant M(d n )0"M(d n ) (or M(d n ~ 1 )0 2 ~M(d n )) electronic 
structure, a wavefunction analysis would enable one to 
decide which model is the most appropriated. On the 
contrary, for intermediate occupation of the oxygen or- 
bital (which is the case in many correlated materials), the 
exact Hamiltonian wave function projection amplitudes 
onto both model spaces are almost the same provided 
that one uses proper orbitals sets. Consequently, the dis- 
crimination between the double exchange and Heisenberg 
model Hamiltonians, which is a matter of numerous dis- 
cussions is simply not possible, neither 
from the spectrum nor from the wave function analysis. 
It is likely that this undecidability will be encountered in 
the interpretation of experimental informations. 
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